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ABSTRACT 

We present a 50 ks observation of the gamma-ray binary LS 1 +61 303 carried out with the ACIS-I array 
aboard the Chandra X-ray Observatory. This is the highest resolution X-ray observation of the source con- 
ducted so far. Possible evidence of an extended structure at a distance between 5" and 12" towards the North 
of LS I +61 303 have been found at a significance level of 3.2a. The asymmetry of the extended emission 
excludes an interpretation in the context of a dust-scattered halo, suggesting an intrinsic nature. On the other 
hand, while the obtained source flux, of Fo.3-10 kev = x 10~ 12 ergs cm" 2 s" 1 , and hydrogen column den- 

sity, Nh = 0.70 ±0.06 x 10 22 cm" 2 , are compatible with previous results, the photon index T = 1.25 ±0.09 is 
the hardest ever found. In light of these new results, we briefly discuss the physics behind the X-ray emission, 
the location of the emitter, and the possible origin of the extended emission ^0.1 pc away from LS 1+61 303. 
Subject headings: stars: emission line, Be — stars: individual: LS I +61 303 — X-rays: binaries — X-rays: 
general — X-rays: individual: LS I +61 303 — X-rays: ISM 



1. INTRODUCTION 

LS I +61 303 is a high mass X-ray binary associated 
with the galactic p l ane va riable radio source GT 0236+610 
( Gregory & Tavlorl 1 1978b which shows periodic non- 
thermal outbursts on average every P or b=2 6.4960±0.0028 d 
dTaylor & Gregory! Il982t iGregoryl 120021) . Optical spec- 
troscopic observations show that the system is composed 
of a rapidly rotating early type BOVe star with a stable 
equatorial decretion disk and mass loss, and a compact 
object wit h a mass between 1 and 4 M w orbiting it every 
~26.5 d dHutchings & Cramptonlll98U ICasares et ail 120051: 
iGrundstrom et alj 120071) . Spectr al line radio observation s 
give a distance of 2.0±0 2 kpc dFrail & Hiellmingl 119911) . 
iMassi et all d200U 120041) reported the discovery of an 
extended jet-like and apparently precessing radio emitting 
structure at angular extensions of 10-50 milliarcseconds. 
Due to the presence of (apparently relativistic) radio emitting 
jets, LS I +61 303 was proposed to be a microquasar. 
However, recent VLBA images obtained durin g a full orbital 
cycle show a rotating elongated morphology (Dh awan et alj 
120061) . which may be consistent with a model based on the 
interaction between the relativistic wind of a non-accretin g 
pulsar and the wind of the stellar companion dDubusll2006l) . 
The MAGIC Cerenkov telescope has recently detected 
LS 1 +61 303 at v ery high energy gamma ray s (>100 GeV; 
Albei-t ltaill2006l). LS I +61 303 and LS 5039 dParedes et alj 
2000; Aharonia n et alj 120051) share the quality of being the 
only two known X-ray binaries that are also GeV emitters. 
Therefore, they can be called gamma-ray binaries. 
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LS I +61 303 has bee n observed w i th di fferent X- 
ray satellites: Einstein dBignami et al.l [l981l). ROSAT 
(Gold oni&Mereghettil 119951: iTavlor et alj 1 19961). ASCA 
dLeahv et al.1 11997b . RXTE dHarrison et all 120001). Bep- 
poSAX, XMM-Newton, and INTEGRAL dSidoli et all 120061: 
IChernvakova et 311 12006). The spectrum could always be fit- 
ted with an absorbed power-law with values in the ranges 
jV H =0.45-0.65 x 10 22 cm" 2 and r= 1.5-1.9. The 2-10 k eV X- 
ray flux shows orbital variability (Par edes et alj|1997h . with 
values ranging from about 5 to 20xl0~ 12 ergs cm" 2 s , the 
maximum occurring at orbital phase ~0.5 (assuming 7o=JD 
2,443,366.775). No spectral lines nor extended X-ray emis- 
sion have ever been reported. Here we present, for the first 
time, Chandra X-ray observations of LS I +61 303, aimed at 
detecting small-scale extended X-ray emission. 

2. CHANDRA OBSERVATIONS AND RESULTS 

We observed LS I +61 303 with Chandra using the standard 
ACIS-I setup in VF mode and 3.24 s time frame during a total 
of 49.9 ks from 2006 April 7 22:30 UT to April 8 12:22 UT 
(MJD 53832.9375-53833.5153). This corresponds to orbital 
phases 0.028-0.050, when a relatively low flux is expected, 
thus minimizing possible pileup effects. 

The Chandra Interactive Analysis of Observations software 
package (CIAO 3.3.0.1) and the CALDB 3.2.2 version have 
been used to perform the reduction of the level 1 event files 
as well as for the spectrum and lightcurve subtraction. The 
reprocessed level 2 event files have been created to account 
for readout effects. The position angle of the readout direction 
of the ACIS-I chip was -62° (positive from North to East). 

We have reduced the data using standard procedures given 
in the Chandra analysis threads. The source region has been 
defined as a 20 pixel (~10") radius circle around the source 
center. For the background, we took a nearby circular re- 
gion of 80 pixel (~40") radius, away from the source and 
the CCD junction. To avoid the presence of possible (very 
faint) sources in the background, the tool celldetect was ap- 
plied. The mkacisrmf and mkarf tools were used to generate 
the RMF and ARF files respectively. The tool dmextract was 
used to generate the source and the background lightcurves 
and spectra. The spectrum was grouped with the tool dm- 



2 



Paredes et al. 



group to 30 counts per energy bin. 

LS 1 +61 303 showed a moderate level of activity, with an 
average count rate of 0.15 counts s" 1 (see Fig.Q~|i, and the data 
was affected by pileup at the ~14% level (see below). The 
background was very low and estimated to be about 1% of 
the source count rate within the source region. The position of 
the X-ray source, determined by using the CIAO package tool 
celldetect, is a = 02 h 40 m 3 1 ?63 ± 0?02, 5 = 6 1 ° 1 3'45 "70 ± 0."3 
(J2000). The errors, as in the rest of this work, are la. The X- 
ray position is consistent with the radio and optical positions. 

2.1. Light curve 

The background subtracted lightcurve of LS 1 +61 303 in 
the 0.3-10 keV band is plotted at the bottom of Fig. Q] Each 
point represents 324 s of data and the error bars represent the 
square root of the total number of counts in each bin divided 
by 324 s. We have also computed the hardness ratio (HR) by 
dividing the count rate in the 1.7-10 keV energy range by the 
corresponding one in the 0.3-1.7 keV range. We show at the 
top of Fig. [T] the HR as a function of time, but with a larger 
bin time of 1024 s to have relatively small error bars. The 
count rate is moderately variable on timescales from several 
minutes to hours, with an average of 0.144±0.029 counts s" 1 
(X 2 ed =l-84). The count rate increases during the observation, 
and a linear fit yields a 30% increase at the end with respect to 
the beginning. Superimposed on this trend is a minifiare 32 ks 
after the start of the observation, when the count rate increases 
by a factor of ^2 over a timescale of ^1000 s, with a total du- 
ration of ^3 ks. This minifiare is accompanied by an increase 
of the HR. This kind of moderate variability was first observed 
with ASCA, which detected so urce variations of a bout a 50% 
on timescales of half an hour dHarrison et al.ll2000l) . Archival 
BeppoSAX observatio ns carried out in 19 97 have also shown 
short-term variability (ISidoli et alJl2006l) . It is interesting to 
note that small amplitude (<10%) radio flares were detected 
in LS I +61 303 wit h a recurrence period of about 5 ks 
dPeracaula et al.l[T997l) . The presence of hour timescale X- 
ray miniflares has also bee n reported in the case of LS 5039 
dBosch-Ramon et ai1l2005l) . A hardness vs. intensity diagram 
for the whole Chandra observation of LS I +61 303 reveals 
a positive trend, the source is harder when it is brighter, but 
with a low significant correlation coefficient of r=0.44. This 
is in contrast to the clear co r relatio n found in previous XMM 
observations by Sido li et all (120061) . 

2.2. X-ray spectrum 

Since the HR is not significantly variable (x 2 ed =1.78), we 
considered the whole data set for the spectral analysis. To 
perform the spectral fits we restricted the energy range be- 
tween 0.3 and 10 keV. A fit to the data with a Bremsstrahlung 
model provided x 2 ed =l-6, and no lines nor a multicolor black- 
body were required to fit the data. An absorbed power-law fit 
yielded x 2 ed =1.14, but a significant excess above 7 keV was 
apparent in the residuals. Due to the count rate and time frame 
used, this is probably due to pile up. The addi tion of the pileup 
model implemented in Sherpa (Davis 2001) to the absorbed 
power law yielded different sets of possible spectral param- 
eter values, all of them equally compatible from a statistical 
point of view. Nevertheless, for some of them the derived 
fluxes were a factor >5 higher than previous measurements at 
similar orbital phases. We found that reasonable results were 
obtained when fixing the event pileup fraction parameter / to 
0.95: Xred = 0-83, grade migration parameter a = 0.33 ±0.09, 
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FIG. 1. — Bottom: Background subtracted lightcurve of LS I +61 303 
obtained with Chandra. The binning time is 324 s and the error bars are 
at the lcr level. Time origin corresponds to 2006 April 7 at 22:33:17 UT 
(MJD 53832.9398), and the observation length is 50 ks. Total count rate 
variations of a factor of two on timescales of a few hours are clearly seen. 
Top: Hardness ratio (HR) between 1.7-10 and 0.3-1.7 keV energy ranges. 
The binning time has been increased to 1024 s to reduce the size of the error 
bars. The HR is not significantly variable and it is poorly correlated with the 
count rate (except during the minifiare). 

N H = 0.70 ±0.06 x 10 22 cm" 2 , Y = 1.25 ±0.09, pileup frac- 
tion 14%, and F . 3 _io ke v = 7-1-H x 10" 12 ergs cm" 2 s" 1 . The 
significance of adding a pileup component to the absorbed 
powerlaw fit was estimated via an F-test: the probability that 
this improvement happens by chance is 0.8%. 

To check the robustness of the results obtained with the 
spectral fitting with the pileup model we have conducted two 
additional tests. First of all we have obtained the spectrum 
of the readout streak, which is not affected by pileup, af- 
ter excluding the 5 pixel diameter region around the source. 
Since it only contains 287 counts, we have fitted the spec- 
trum using low count number statistical weights. Fixing A^h 
to 0.7 x 10 22 cm" 2 we have obtained T = 1.21 ±0.14, which 
is consistent within errors to the value given above. An addi- 
tional test has been conducted by fitting the spectrum of the 
source wings, which are not affected by pileup. We have built 
a spectrum with data beyond a radius of 1.5" from the source 
center (i.e., excluding the inner 3x3 pixels, those most af- 
fected by pileup) and fitted it with and without fixing A^h, ob- 
taining in both cases a photon index harder than the one ob- 
tained with the pileup model. A similar test with data beyond 
2.0" has provided similar results. These tests clearly suggest 
that the photon index value obtained with the pileup model is 
reliable, and that the source was intrinsically hard during this 
observation. Moreover, the estimate of the total source flux 
using the fitted flux of the source wings and the fraction of 
the simulated PSF (see next section) counts contained in the 
wings is compatible for /=0.95. This is not the case for the 
lower values of / that implied the much higher source fluxes 
discussed in the previous paragraph. These results give us 
confidence on the results obtained with the pileup model. 

Regarding the emission of LS I +61 303 itself, the flux 
we report here, Fo.3-10 kev = 7.1*^| x 10~ 12 ergs cm" 2 s" 1 , is 
compatible within errors to that obt ained during a short XMM 
observation also around phase (IChernyakova et all 120061) . 
However, the photon index Y = 1 .25 ±0.09 found in the Chan- 
dra data is the hardest ever found, only compatible with any 
previous value at the 3er level and harder, at a ^5cr level, than 
the XMM one at the same phase (Y = 1 .78 ± 0.04). It is out- 
side the scope of this work to put forward a physical model 
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FIG. 2.— Le/f: Image of LS I +61 303 obtained with Chandra in the 0.3-10 keV energy range. The image size is 24"5 X 24"5 (the pixel size is 0"49 X 0"49). 
A Gaussian kernel of three pixels in radius has been used to smooth the image. The intensity scales logarithmically. The contours correspond to seven logarithmic 
steps in the range 0.2-300 counts pixel -1 . Extended X-ray emission is probably present at a distance of 5-12" from the center. Middle: Same as before but 
for the PSF. Right: Same as before but for the residuals obtained after subtracting the PSF image to the source one. The intensity scales here via a histogram 
equalization. The extended emission at distances beyond 5", in the form of several spots towards the North-East, is apparent. The central residuals at ~2" radius 
occur within the core of the PSF, and no conclusion can be drawn because of pileup. 



for the spectral state of the source. We will just point out 
that the X-ray photon index can be harder than the canoni- 
cal value of 1.5, which corresponds to a electron population 
with power-law index of 2. This can occur for different rea- 
sons, for examp le a Klein-Nish ina dominated steady particle 
population (see Derishev 2006] for an extended discussion). 
A very hard spectral component (e.g. inverse Compton or 
Bremsstrahlung), produced by the lower energy part of the 
particle spectrum with a high low-energy cutoff, may also 
harden the synchrotron spectrum. The lack of simultaneous 
data at higher energies does not allow us to distinguish be- 
tween the two hypotheses. 

We finally selected the data corresponding to the minifiare 
and obtained a spectrum for it, which had poor statistics and 
was especially affected by pileup. A fit to this spectrum 
yielded a photon index compatible within (large) errors to that 
of the whole observation spectrum. 

2.3. Imaging 

A preliminary visual inspection of the source image re- 
veals the possible existence of extended emission. To better 
study this issue we have used ChaRT and MARX to simulate 
the point spread function (PSF) accounting for pileup effects. 
This PSF has the same counts as the source, which are dis- 
tributed in energy according to the source spectrum including 
pileup, being spatially distributed in the image considering 
pileup as well 6 . We show in Fig.[2]the images corresponding 
to the source, the PSF, and the residuals obtained after sub- 
tracting the PSF from the source, all of them smoothed with a 
Gaussian kernel of three pixels in radius to enhance the faint 
emission. As can be seen in Fig. |2}left, there appears to be 
extended X-ray emission around LS I +61 303 at distances 
between 5" and 12" from the source center. The subtrac- 
tion of the nearly circular PSF, shown in Fig. |2]-middle, from 
the source image yields the residuals shown in Fig. |2}right. 
There seems to be a region with an excess of counts at several 
arcseconds from the image core and nearly half-surrounding 
it. This excess region, already apparent in the source image, 

6 We note that, while the pileup model included in MARX is valid for qual- 
itative predictions of the effects of pileup on the PSF, it has not been verified 
for image reconstruction (section 7.3 of the MARX 4.0 Technical Manual). 
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FIG. 3. — Surface brightness distribution of LS I +61 303 (red filled circles) 
compared to the PSF (dashed blue line). The excess above ~5" is apparent 
in this radial profile, even if the extended emission is asymmetric. 

presents some spots and extends mainly towards the North- 
East. 

To check the reliability of the 5-12" excess we have cre- 
ated radial profiles of the surface brightness for the source 
and the simulated PSF with the tools dmextract and dmtcalc. 
We show the radial profiles with bins of 1" in Fig. [3] where 
hints of excess are seen in the region 5-12" away from the 
center. Moreover, we note that a comparison with the un- 
scattered point-like blazar source PKS 2155-304 located at 
high-galactic latitude (following the procedure described in 
iGallo & Fended l2002l) provided an even clearer departure of 
LS 1+61 303 from a point-like source. Nevertheless, we pre- 
fer to show here the most conservative results obtained with 
the PSF, simulated at the chip position and with the observed 
spectrum of LS 1 +61 303 including pileup effects. Next we 
computed the excess of source counts over the PSF in an an- 
nulus with inner and outer radii of 5" and 12. "5, and obtained 
58±18 counts (3.2<r, or single trial probability of 1.4 x 10" 3 ). 
We have looked for excesses in a circle with radius 20" cen- 
tered on the source, without considering the central 5" radius 
circle, thus yielding an area of 7r(20 2 -5 2 )" 2 . Dividing this by 
the area of the annulus where we obtained the 3.2cr excess we 
find that the number of trials is 2.9, yielding a probability of 
4. 1 x 10~ 3 , which corresponds to 2.9a. However, since the ex- 
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tended emission appears to be asymmetric, we have computed 
the excess dividing the 5-12"5 annulus in four quadrants. Al- 
though there is no excess towards the West and South direc- 
tions, the excess is of about 1.5c to the East, and reaches 3.8<r 
to the North (or single trial probability of 1 .4 x 10~ 4 ). We note 
that this excess is not located in the path of the readout streak. 
The number of trials is now 1 1.6, providing a post-trial prob- 
ability of 1.6 x 10~ 3 or a 3.2a detection. The excess in the 
Northern quadrant of the 5-12"5 annulus, with a projected 
distance of 0.05-0.12 pc from the center of LS I +61 303, 
represents 0.5% of the source counts, implying (with T=1.25) 
a flux of ^4x 10~ 14 ergs cm" 2 s" 1 and Lx ^2x 10 31 ergs s" 1 . 

3. DISCUSSION 

The Chandra observations reported here show that 
LS I +61 303 presented moderate fluxes of -Fb 3-10 kev = 
7.l!i jxl0" 12 ergs cm 2 s 1 , as expected around orbital phase 
~ (ISidoli etalJ (| 200 6Y). The source appeared to be moder- 
ately variable on timescales of ^1000 s. The hydrogen col- 
umn density was A^h = 0.70 ± 0.06 x 10 22 cm" 2 , similar to the 
values found in the past. The photon index was 1.25 ±0.09, 
harder than the previously measured values in the range 1 .5- 
1.9. Finally, an excess between 5 and 12"5 towards the North 
of LS 1+61 303 has been detected at the 3.2er level (post-trial). 

Chandra observations of X-ray binaries have provided the 
most detailed X-ray images of these sources. Extended circu- 
lar emission, produced by ISM dust grai ns along the line of 
sight that scatter mainly the soft X-rays ( Prede hl & Schmitl 
1995), has been found in most cases (IXiang et alJ 120051) . 
However, the possible extended structure around LS I +6 1 303 
we report here is asymmetric, suggesting its intrinsic nature. 
Plausible explanations of this diffuse radiation could be ther- 
mal Bremsstrahlung produced by hot and dense ISM close 
to the source, perhaps heated and compressed by a large-scale 
outflow originated in LS I +6 1 303, and synchrotron or inverse 
Compton emission produced by non-thermal particles accel- 
erated at the termination region of the aforementioned out- 
flow. The first possibility would require a medium with den- 
sities significantly higher than those inferred from the ISM ab- 
sorption towards LS I +61 303. If non-thermal, the extended 
emission would more likely be of synchrotron origin due to 
its shorter timescales to radiate at a few keV than those of 



inverse Compton or relativistic Bremsstrahlung. In any case, 
the asymmetric distribution of counts could be related to the 
main direction of the outflow and/or to ISM inhomogeneities. 
Regardless of the morphology, we can say that efficient en- 
ergy transport is probably taking place up to distances of at 
least ~0.1 pc in LS 1 +61 303. However, only deeper high- 
resolution observations can help to better detect the excess 
(currently at a 3.8tr level, 3.2c post-trial) and determine the 
precise morphology and underlying physics. 

A final comment on the hydrogen column density appears 
necessary. The A^h value inferred from spectral fits to all avail- 
able X-ray observations is in the range 0.45-0.70 x 10 22 cm -2 , 
while the A^h value of the ISM inferred from UV/optical 
absorption is in the range 0.45-0.60 x 10 22 cm" 2 . Therefore, 
the detected X-rays suffer small intrinsic absorption, if any. 
Moreover no significant A^h orbital vari ations have ever been 
detected (as pointed out previously by Leahv et al. 1997] after 
ASCA observations). In conclusion, the apparent lack of 
significant soft X-ray absorption and the constancy of A^h 
along the orbit are difficult to be explained considering the 
slow and dense material expelled from the companion star. 
The X-ray emitter should be placed away from the binary 
system to produce such observational results, in contrast 
with the expectations from classical colliding wind and 
microquasar scenarios. This issue certainly requires further 
investigation. 
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